Revisiting Langer-Ambegaokar-McCumber-Halperin theory of resistive transitions in 
one-dimensional superconductors with exact solutions. 
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We present an important correction to the Langer-Ambegaokar-McCumber-Halperin theory for 
the resistive state of a ID superconductor. We establish that the identification of the saddle on the 
free energy surface over which Langer and Ambegaokar had claimed the system to move in order to 
form thermally excited phase slip centres is wrong. With the help of an exact solution we show that 
the system has to overcome a similar free energy barrier but can actually have vanishing amplitude 
of superconducting phase at a point unlike the Langer-Ambegaokar solution. 
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In a ID superconductor, kept below the critical tem- 
perature Tq, a current driven transition to a normal (N) 
phase from the superconducting (SC) phase has been 
seen to occur in experiments as early as in 1967 [l|, 
Such transitions show a finite width in temperature and 
current strength within which SC and N phases coex- 
ist. A theory was proposed by Langer and Ambegaokar 
(LA) which was subsequently improved by McCum- 
ber and Halperin (MH) jj], to explain thermal fluctua- 
tions induced transitions through metastable states. Th 
Langer-Ambegaokar-McCumber-Halperin (LAMH) the- 
ory traces the origin of localized N phases in supercon- 
ducting ID (width is smaller than the coherence length £) 
samples in the formation of phase slip centres (PSC) - an 
idea which was originally put forward by Little 5] . Over 
last 40 years, the LAMH theory successfully accounts for 
the resistivity vs temperature plots observed in this (so- 
called) resistive regime except for some stretches at the 
lower temperature end of the resistivity vs temperature 
plot. At this end (near the Tc), where the resistivity 
vanishes and the system moves into the SC phase, some 
mismatch with the theoretically predicted values have 
been observed in the early experiments which were be- 
lieved to be effects of imperfections at the contacts at the 
ends of the sample @ . These deviations have also been 
seen to be system specific particularly lending support to 
relating their origin to contact imperfections. Classic ex- 
periments done by Webb-Warburton 0] and Newbower et 
al [1] showed the excellent applicability of LAMH theory 
to experimental results. Since then, the LAMH theory 
has remained the basic tool for dealing with the resis- 
tive regime of the ID superconductor near the Tc regime 
where quantum effects are negligible. 

Such a successful theory, however, has a particular la- 
cuna which is a glaring inconsistency in it. The SC order 
parameter is a complex number and in ID it looks like 
a spiral wound around the wire. The (LAMH) theory is 
based on a suggestion by Little, that, the amplitude of 
the SC phase has to locally vanish at a point along the 
length of the ID sample in order to have a turn added or 
removed from it 51. Based on this criterion, one looks for 



an amplitude modulated solution of the SC phase where 
the amplitude vanishes at least locally (PSC) allowing the 
SC order parameter to add or remove turns (change its 
wave number q). Thus, the total phase along the length 
of the sample can change. An applied voltage across the 
length of the sample, on one hand, keeps adding turns 
to the SC order parameter (increases wavenumber) to 
increase the phase. On the other hand, thermal fluctua- 
tions make the system access unstable amplitude modu- 
lated states like PSCs to give up phase in the middle of 
the sample. The LA theory basically rests upon a balance 
between these two. Although, this vanishing of the am- 
plitude is crucial for such a phenomenon, the LA solution 
which corresponds to the passage of the system through 
a free energy saddle never goes to zero anywhere. Inter- 
esting to note that, the LA solution can never actually 
go to zero because of infinite free energy cost it has to 
incur according to the LA calculations of the free energy 
and thus it is absolutely against the demand of Little's 
criterion. LA theory simply claims that by some other 
fluctuations the solution that comes close to zero would 
become zero which is wrong in view of the free energy 
surface as described by LA. Despite this glaring mistake, 
the LA theory works fine because in reality (on correct 
free energy surface) the free energy barriers as calculated 
by LA and the true one (which will be shown) arc the 
same. 

In the present paper, we will first clearly identify the 
error in the identification of the saddle on the free energy 
surface in LA theory and demonstrate why the good old 
expression for the variation of amplitude at the bundary 
of a superconductor under no-field condition holds as an 
exact solution in the case of LA model as well. Particu- 
larly on the basis of this exact solution we will also cal- 
culate the corresponding chemical potential profile and 
the current profile through the sample. We would show 
that the divergence of the current at the point where 
the amplitude vanishesh goes against the LA demand of 
other fluctuations taking a close to zero amplitude ac- 
tually to zero because, the LA effective free energy can 
not allow that, we would also argue that the addition 
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or subtraction of turns in the middle of the supercon- 
ducting sample is not a process that takes place just at 
the point the amplitude vanishesh, but, should take place 
rather continuously through the process of formation and 
relaxation of PSCs. 

The dynamics of a ID superconductor in the resis- 
tive regime is given by time dependent Ginzburg-Landau 
(TDGL) equation [i| as 



(tpt + ifiip) = tpxx + (a- ^\ip\ 2 )ip 
j = Im{ip* v tp) - Vx- 



(1) 
(2) 



In this model, ip is the superconducting order parameter 
which is complex valued. The system being ID, x is 
distance along the wire from some arbitrary origin. Let 
us consider, the length of the wire is L with cross sectional 
area a. The [i is electrochemical potential which can be 
considered as the order parameter of the N phase within 
the scope of G-L phenomenology. The constants a and (3 
measures the free energy density differences of the SC and 
N phase as g n — g s = a 2 /2(3 = a(AT 2 ) where a is another 
constant that depends upon density of states at Fermi 
surface iV(0) and Boltzmann constant as a — 4.7N(0)k%. 
The j is current (density) through the sample and the 
suffix t and x of ip and /i indicate of partial derivatives. 
Eq.l has two stationary solutions. (1) if> = 0, /i 
.;•;/ (for constant j) which is the normal state and (2) 
ip = Ae iqx , A 2 = (a - q 2 )/(3, j = A 2 q, which is the 
superconducting state when fj, = (3j|. The SC order 
parameter can be visualized as a spiral wound around the 
x-axis (along the wire). If there are N turns along the 
length L, there is a total phase difference <\> = 2ttN along 
L. Thus, the wave number q of the SC phase is a measure 
of number of turns the system has on a given length L, 
since, q — 2nN/L. The expression for the corresponding 
G-L free energy, when the steady state solution is purely 
superconducting, is given by 



F = L(i[(q 2 - a)A 2 



A 4 



(3) 



where a is the cross section of the wire. The form of F 
clearly indicates that, the SC states with smaller q val- 
ues are energetically favoured. In other words, the spi- 
ralling SC order parameter would tend to lose its turns 
to go to a lower free energy state. Putting the ansatz 
■ijj = A(x,t)e l ^ x>> where d<f>(x)/dx = q(x) in Eq.l and 
separating the real and imaginary parts (where we con- 
sider A(x, t) and <p(x, t) real functions) we get 

f) A f) 2 A 

^-^-(a^q 2 )A + ^=0 (4) 
A dq dA A 

along with these two equations mentioned above one has 
to take into account the Eq.2 to get a complete picture 
of the affairs. 



Before we go into the exact solution and calculation 
of the barriers, let us first have a look at what had gone 
wrong in the LA calculations of the barrier. The station- 
ary Eq.4 can be rewritten in the form 



d 2 A 

dx 2 



(«-g 2 ) 



A 2 
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The above expression clearly shows that the effective po- 
tential U goes to zero at A = 0. Whereas, in the LA 
theory, because of replacing the q 2 A 2 /2 with — j 2 /2A 2 
in the expression of the U, it diverges at A — for all 
nonzero j. This replacement of q 2 A 2 /2 with j 2 /2A 2 is 
clearly wrong because the amplitude in the above ex- 
pression is not the super-current amplitude rather its a 
modulated form of it and the current density in such a 
case is not the super-current density corresponding to the 
constant amplitude SC phase. 

Fig. la and b show a schematic comparison of the 
U as taken by LA and us. In our case, the A = 
point would be visited by the trajectory from the point 
A 2 = (a — q 2 )//3 where the state is metastable at the cost 
of an increase of the velocity equivalent 4^ . Interesting 
to note that, the point A 2 — (a — q 2 )/ft is the peak of the 
U similar to the one considered in LA theory. But, the 
U of LA does not allow the system to reach the A = 
because that would require an infinite storage of kinetic 
energy equivalent ( j^r) 2 /2. The actual Ginzburg-Landau 
free energy containing a gradient square term would then 
have to keep an infinite free energy storage somewhere in 
the form of a diverging amplitude gradient in order to 
allow the system to ever reach the A = point. But, the 
constant amplitude SC phase does not allow for such a 
storage. To circumvent this inconsistency LA proposes 
this barrier to be a saddle such that the system can es- 
cape through other dimensions of the space in which the 
free energy is defined. A replacement of q by j and A still 
keeps the space two dimensional and there we see that the 
j has to fall faster than A and necessarily vanish to keep 
the free energy divergence free as the A vanishes. This 
goes against the very concept of current driven origin of 
the resistive state. Since, the LA calculations mean that 
the A = is never accessible for a nonzero j, the subse- 
quent claim of LA theory that other fluctuations would 
make LA solution which goes closer to zero actually reach 
zero is untenable. As a consequence of a wrong replace- 
ment of the wave number the saddle identified is not the 
lowest barrier and that would be clear in the following 
where we would show a smaller barrier. 

Note that, at a constant q the Eq.5 gets a form which 
had been arrived at by one of the authors @, [Io| pre- 
viously by a separation of length scales. Eq.5 gives the 
corresponding /i profile to an amplitude modulation ob- 
tained from Eq.4. In what follows we will stick to this 
constant q scenario because, a. Eq.4 that actually ad- 
mits the amplitude modulations is independent of any 
variation of q, b. even with constant q we will be able to 
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(a) 




(b) 



FIG. 1. schematic diagrams for the comparison of ef- 
fective potentials on arbitrary scales: (a) U as consid- 
ered by LA and (b) U that we consider here. 



show an exact form of the PSC which has a smaller free 
energy barrier than that proposed by LA and that serves 
our present purpose and c. this is a good approximation 
because, normally a single turn gets added or removed 
at the formation of the PSC and that changes the q by 
an order 1/L which is quite small. Considering the time 
independence of the amplitude, Eq.4 admits a solution 
A = Aq tanhx/v^ with the couple of conditions - (1) 
(a - q 2 ) - l/£ 2 = and (2) (3A 2 - l/<f = where £ 
is the bulk GL coherence length. These conditions im- 
mediately identify that Aq 2 = (a — q 2 ) / (3 which is the 
standard amplitude of the steady superconducting state. 
From the above conditons we recover the standard bulk 
relationship between the a and £ also. The above men- 
tioned solution is a well known one which one observes 
at the boundary a bulk superconductor in no-field condi- 
tions. The LA free energy denies it simply because of the 
divergence of the effective free energy at zero amplitude. 

The corresponding /i profile of the superposed normal 
phase comes out to be /i = (2q/y/2£)(l/ tanhx/%/2^ — 
tanhcc/v^O- This is the shouted after exact solution 
which would corresponds to an energy barrier same as 
that in LA theory. The barrier height would become even 
smaller with a nonzero q as is expected from GL free en- 
ergy. Corresponding to our exact solutions for A and fi 
the current can be evaluated from Eq.2 for a constant 
q. Its important to note that, this current has a value 



A\q far from the origin and it diverges at the origin as 
1/ tanh 2 x/\/2~£- So, it has a part which is grows as 1/A 2 
where A — > and there is no question of the current 
falling faster than A for the divergence in the effective 
free energy U of LA to remain analytic. So, according to 
the free energy expression of the LA this exact solution 
of the system does not qualify for overcoming the free 
energy barrier but so far done experimental studies basi- 
cally confirm having such solutions that vanish locally. 

At the limit the super-current density j — ¥ the 
very complicated looking LA solution through the LA 
mentioned saddle of the free energy surface actually 
overcomes a free energy barrier somewhat bigger than 
AF = (8V2/3K B T)(g n - g s )a£ @, 1, Ej. The free 
energy barrier corresponding to our exact PSC solu- 
tion A — Aq tanh x/y2^ can be calculated by putting 
ijj = Aq tanh x/y2^ e lqx and A$e %qx respectively in the ex- 
pression of the free energy (corresponding to Eq.l) shown 
below and subtracting the latter free energy from the for- 
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dx[\ v^l 2 -c*M 2 



(7) 



If we make use of the conditions (1) and (2) as mentioned 
above in connection with the derivation of the PSC we 
can immediately show that the free energy barrier is ex- 
actly calculated as 



A F = (8V2/3K B T)(g n - 5 >£ - 



26y/2a 2 28ct£ 4 
(8) 

We can see that the q = that corresponds well with the 
j = A 2 q = condition, the free energy barrier is clearly 
the same as the LA value when the relevant length scale 
of the amplitude modulation is equal to the bulk GL co- 
herence length £. Taking roughly q l/£ one recovers 
the AT 2 / 3 power law in all the terms of the above men- 
tioned free energy expression. Note that, consideration 
of a nonzero q even reduces the free energy barrier. Im- 
portant to note that q — 0, actually sets fi = which is 
a violation of the relation 



Aire , d A arqib 
A7= 

h dt ' 



(9) 



where AV is the applied voltage drop across the system. 
In literature, the limit at which the AT becomes ATc* 
is 3 1 3c ~ > where jc — 2a 3 / 2 /3 v / 3/3- This limit can be 
seen as q — > limit also, but, q cannot actually vanish 
for the theory to be consistent and so is j. Taking into 
account that, a nonzero q is essential one can easily infer 
that a nonzero A and j, however small, should accom- 
pany a PSC and that the LA solution of the PSC would 
actually encounter an enormous barrier. 

Let us have a discussion on the implication of the 
present results. The LAMH theory is the most useful 
and accepted framework in understanding classical PSC 
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induced origin of the resistive regime of the ID super- 
conductor. It matches pretty well the experimentally 
obtained resistivity vs temperature plots. Our present 
analysis helps solve the riddle as to how such a success- 
ful theory can actually not predict a proper PSC solution 
and establishes LAMH theory on even stronger founda- 
tion after more than 40 years of its introduction. We 
clearly show that, it was the wrong identification of the 
free energy barrier which resulted in this discrepancy in 
the value of the ATq by finding out an exact solution of 
the same dynamical equation used by LAMH and corre- 
sponding free energy barrier. The LA theory, however, 
uses the correct numerical value despite identifying the 
wrong barrier because it had stuck to the conservation of 
energy in getting the saddle. Here, we actually show that 
no such saddle exists because the current density can not 
go to zero at the point as the amplitude does. We have 
shown that the good old solution for the amplitude mod- 
ulations at the boundary of a bulk superconductor holds 
perfectly good in the LA case as well. 

Our present analysis apart from resolving this PSC 
riddle also clarifies and sheds light to a few other points. 
First of all, we have considered here a constant q during 
the PSC formation. This not only simplifies the calcula- 
tions but comes out to be correct because, in general only 
a turn addition or removal happens by the formation of a 
PSC. Moreover, as we have already shown that, there is 
no reason to think of a coupling between the amplitude 
and the wave number in the modulated amplitude phase 
just as in the SC phase and that has severe consequences. 
It has to be explored what other forms of PSC like or 
other solutions exist with a spatially varying q. Never- 
theless, in the present context, with a constant q we are 
able to point out the discrepancy in the calculations of 
the saddle in LA theory. Our solution which satisfies the 
Little's criterion of amplitude vanishing is clearly having 
the same barrier to overcome as LA's saddle and in that 
sense is a better representative for PSC. Note that, our 
exact solution is unstable for tanh (x/\/2£) < 1/3 to 
any infinitesimal uniform perturbation. So, fluctuations 
would start growing at the core of the PSC and would 
make it relax back to the constant amplitude SC phase. 

Another important point to note is that, the claim that 
the LA solution approximates to a form proportional to 
tanh|a;| ll| is somewhat in conflict with the Eq.4 be- 



interesting to note that, the Eq.l is invariant under a con- 
stant global phase shift. As a result, such a phase shift 
should not be energetically unfavourable and can happen 
with a locally vanishing amplitude particularly because 
the single valued nature of the SC order parameter is not 
compromised. But, it marks a little departure from the 
Little's proposition that a phase change of 2ir or its mul- 
tiple would actually happen right at the vanishing of the 
amplitude. Rather, we see here that a phase change of 
half of that required to add/lose a turn is happening at 
the time of the formation of the PSC and the other half 
has to happen during the relaxation of the PSC to keep 
the order parameter single valued. One way of looking 
at this scenario would be that, as the PSC forms, it pro- 
duces strain locally to the turns and that strain relaxes 
by rotating one half of the spiral. By the time the PSC 
has occurred, half of the spiral is rotated by an angle 
7r. Now a relaxation of the PSC might find it energet- 
ically favourable or retaines some memory to continue 
effectively rotating the same half in the same direction 
so that when the relaxation of the PSC is complete the 
system loses/gains a turn. A detailed investigation using 
the dynamics could be revealing. We conclude by saying 
that, we have exactly solved the TDGL model for the 
resistive state of an one dimensional superconductor to 
explicitly show a similar energy barrier as the LA saddle. 
On the basis of that we identify the error in the LA cal- 
culation of energy barrier. We put the LAMH theory on 
a even firmer basis. 



cause, this form would produce a delta function contri- 
bution at the origin in Eq.4 resulting in an invalidity of 
such a form as a solution. This is why our exact solution 
is of the form tanh (x/ which induces a global phase 
shift of 7r on one half of it unlike the tanh Id form. Its 
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